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[l] The 8.2 ka event was the last deglacial abrupt climate 
event. A reduction in the Atlantic meridional overturning 
circulation (AMOC) attributed to the drainage of glacial 
Lake Agassiz may have caused the event, but the freshwater 
signature of Lake Agassiz discharge has yet to be identified 
in 6 18 0 of foraminiferal calcite records from the Labrador 
Sea, calling into question the connection between freshwater 
discharge to the North Atlantic and AMOC strength. Using 
Mg/Ca-paleothermometry, we demonstrate that ~3°C of 
near-surface ocean cooling masked an ~1.0%o decrease in 
western Labrador Sea (WO of seawater concurrent with 
Lake Agassiz drainage. Comparison with North Atlantic 
6 ls O of seawater records shows that the freshwater 
discharge was transported to regions of deep-water formation 
where it could perturb AMOC and force the 8.2 ka event. 
Citation: Hoffman, J. S., A. E. Carlson, K. Winsor, G. P. 
Klinkhammer, A. N. LeGrande, J. T. Andrews, and J. C. Strasser 
(2012), Linking the 8.2 ka event and its freshwater forcing in the 
Labrador Sea, Geophvs. Res. Lett., 39, L18703, doi:10.1029/ 
2012GL053047. 

1. Introduction 

[ 2 ] The 8.2 ka event is recorded in numerous paleoclimate 
records throughout the Northern Hemisphere [Alley et al, 
1997; Alley and Agustdottir, 2005; Rohling and Pdlike, 
2005]. Proxies of Atlantic meridional overturning circula- 
tion (AMOC) (e.g., sortable silt, benthic S 13 C) indicate a 
reduction in AMOC during the 8.2 ka event, suggesting the 
involvement of AMOC variability in this abrupt climate 
event [Ellison et al., 2006; Kleiven et al., 2008; Hoogakker 
et al., 2011]. Earlier abrupt climate events with a similar 
climatic footprint are linked to increased freshwater dis- 
charge to the North Atlantic and subsequent reductions in 
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AMOC strength [Clark et al., 2001], leading to the hypoth- 
esis that increased North Atlantic freshwater discharge like- 
wise caused the 8.2 ka event [Klitgaard-Kristensen et al., 
1998; Barber et al., 1999; Clark et al, 2001; Ellison et al., 
2006; LeGrande et al., 2006; Came et al., 2007]. 

[ 3 ] During the last deglaciation, large proglacial lakes 
formed between the retreating Laurentide Ice Sheet margin 
and the isostatically depressed landscape [Teller etal, 20021. 
Lake volume reached a maximum (40,000-151,000 km 3 ) 
when glacial Lake Ojibway combined with Lake Agassiz just 
prior to ^8.5 ka [Teller et al., 2002; Clarke et al., 2004]. 
Upon collapse of the Laurentide Ice Sheet over Hudson Bay 
shortly after ~8.5 ka [Barber et al., 1999], Lake Agassiz 
drained through Hudson Bay into the Labrador Sea at an 
estimated discharge of ~5 Sverdrups (Sv, 10 6 m 3 s _1 ) in one 
to several ^6-month -long flood(s) [Andrews et al, 1995, 
1999; Kerwin, 1996; Clarke etal., 2004; Ellison et al., 2006; 
Hillaire-Marcel et al, 2007; Lajeunesse and St-Onge, 2008; 
Clarke et al., 2009; Roy et al., 2011; Lewis et al., 2012], 
In addition to the freshwater forcing from the flood(s), col- 
lapse of the Laurentide Ice Sheet also routed Lake Agassiz 
runoff through Hudson Strait [Clark et al., 2001; Hillaire- 
Marcel et al., 2007; Clarke et al., 2009] for several hun- 
dred years at a discharge of 0.13 ± 0.03 Sv [Carlson et al., 
2009]. If this freshwater reached regions of deep-water for- 
mation in the North Atlantic, it could have capped surface 
water and disrupted the AMOC [Ellison et al., 2006; Kleiven 
et al, 2008; Hoogakker et al., 2011], The relatively close 
timing of the 8.2 ka event [Alley et al., 1997; Alley and 
Agustdottir, 2005; Rohling and Pdlike, 2005] and the drain- 
age of Lake Agassiz suggests a causal relationship where 
lake drainage and attendant runoff routing could have forced 
this abrupt climate event [Klitgaard-Kristensen et al., 1998; 
Barber et al., 1999; Clark et al., 2001; Teller et al., 2002; 
LeGrande et al., 2006; Hillaire-Marcel et al., 2007; Clarke 
et al., 2009; Carlson et al, 2009; Lewis et al, 2012], 

[ 4 ] Sediment and faunal records from James and Hudson 
Bays, Hudson Strait, and the Labrador Sea support the 
drainage of Lake Agassiz around 8.2 ka [Andrews et al., 
1995, 1999; Kerwin, 1996; Barber et al, 1999; Hillaire- 
Marcel et al., 2007; Lajeunesse and St-Onge, 2008; Roy 
et al, 2011; Lewis et al, 2012], and Hudson Strait geo- 
chemical records document the routing of Lake Agassiz 
drainage basin runoff during the 8.2 event [Carlson et al., 
2009], However, no direct evidence of a freshwater signal 
has been found in Labrador Sea d 18 0 of foraminiferal calcite 
(6 18 O c ) coeval with the drainage of Lake Agassiz [Keigwin 
et al., 2005; Hillaire-Marcel et al., 2007, 2008]. One expla- 
nation for the lack of a d 18 O c signal is that the lake drainage 
was too short lived to be recorded by planktonic foraminifera 
[Andrews etal., 1999; Hillaire-Marcel etal., 2007], although 
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Figure 1. Locations of cores HU87033-017 (54.6°N, 
56.2°W; this study), MD99-2203 (35.0°N, 75.2°W) 
[Cleroux et al., 2012], MD99-2251 (57.4°N, 27.9°W) 
[Ellison et al., 2006], ODP Site 984 (61.0°N, 25.0°W) 
[Came et al., 2007], and RAPiD-12-lK (62.1°N, 17.8°W) 
[Thornallev et al., 2009] indicated by circles with color 
showing observed b 18 O sw anomaly during the 8.2 ka event. 
Black circles indicate locations of decreases in <5 18 0 of cal- 
cite [Keigwin et al., 2005]; <5 l8 O sw changes for these records 
are not displayed because temperature effects have not been 
assessed. Also shown is the GISS ModelE-R simulated 
change in surface ocean b 18 O sw one decade after the addition 
of 5 Sv of freshwater in 0.5 years [LeGrande and Schmidt, 
2008] (<5 18 0 of — 25%o) [Hillaire-Marcel et al., 2008] into 
Hudson Strait with an already weakened convection in the 
Labrador Sea (run WK 5 Sv x 0.5 year) [LeGrande and 
Schmidt, 2008]. Note that the foraminifera records reflect 
b 18 O sw changes at the preferred depth of foraminifera, which 
may be deeper than the modeled surface ocean b 18 O sw 
changes. 

a longer b 1 8 O c decrease should be evident from the breakup 
of the Laurentide Ice Sheet over Hudson Bay [Andrews etal, 
1995, 1999] and the addition of the Lake Agassiz drainage 
basin to Hudson Bay [Clark et al., 2001; Hillaire-Marcel 
et al., 2007; Carlson et al, 2009]. Another explanation is 
that the discharged freshwater may have had little impact on 
Labrador Sea <5 lx O c because Lakes Agassiz and Ojibway 
<5 18 0 of approximately — 25%o was only marginally more 
negative than present-day runoff to James and Hudson Bays 
of — 20%o [Hillaire-Marcel et al., 2008]. However, even 
without a change in runoff b 18 0, an increase in freshwater 
discharge as expected during Lake Agassiz drainage would 
still cause a decrease in Labrador Sea b 18 0 of seawater 
(b 18 O sw ) due to the increased volume of 18 0-depleted fresh- 
water [Barber et al., 1999; LeGrande et al., 2006]. Attrib- 
uting the cause of the reduction in AMOC during the 8.2 ka 
event to the drainage of Lake Agassiz has also proven diffi- 
cult, because Lake Agassiz discharge may have been mostly 
trapped against the Labrador shelf in a buoyant surface flow 
of the Labrador Current and not transported to regions of 
deep convection in the open ocean [Wunsch, 2010; Condron 
and Winsor, 2011]. These outstanding issues ultimately 
imply that freshwater discharge may not have caused the 
8.2 ka event, raising the possibility that alternative 


mechanisms, like reduced solar activity, forced this and other 
abrupt climate events [Alley and Agustdottir, 2005; Rohling 
and Palike, 2005]. Because b 18 O c is sensitive to changes in 
surface temperatures and local <5 18 O sw [e.g., LeGrande et al., 

2006] , we hypothesize that the magnitude of near-surface 
cooling during the 8.2 ka event masked a b 1 8 O sw decrease in 
Labrador Sea b 1 x O c from the drainage of Lake Agassiz. 

2. Methods 

[5] We test our hypothesis with core HU87033-017 from 
the Cartwright Saddle (54.62°N, 56.18°W, 514 m water 
depth) (Figure 1 and auxiliary material, Text SI, section SI), 
which is located in the pathway Lake Agassiz discharge 
would have taken upon exiting Hudson Strait [Condron and 
Winsor, 201 1]. 1 We investigate early Holocene sediment 
discharge from the Laurentide Ice Sheet with changes in 
weight percent coarse fraction (>63 pm). We document 
western Labrador Sea calcification temperatures (CT) of 
Neogloboquadrina pachyderma (sinistral) with flow-through 
Mg/Ca-paleothermometry (Text SI, section S4) [Klinkhammer 
et al., 2004], N. pachyderma (s) is a pycnocline dwelling 
subpolar to polar planktonic foraminifer (Figure SI and 
Text SI, sections S2 and S7). Combining ourMg/Ca CT record 
with the existing b 18 O c record of Andrews et al. [1999] from 
HU87033-017, we calculate b 18 O sw (Text SI, sections S4-S8). 
Due to the low abundance of planktonic foraminifera tests, the 
core chronology is based on eight benthic, reservoir-corrected 
(450 years), calibrated 14 C dates of Andrews et al. [1999] 
(Table SI and Text SI, section S3). 

3. Results 

[6] Our new grain size record for this core identifies 
an abrupt ~20% increase in weight percent coarse fraction 
(>63 //in) dated at ~8.5— 8.3 ka (Figure 2b), similar to the 
~10% increase in the existing >125 pm record of Andrews 
et al. [1999]. This layer lies directly above the final degla- 
cial increase in weight percent calcite and dolomite 
[Andrews et al., 1999] (Figure 2b). The <63 pm fraction 
of these layers has a reddish hue, indicative of sediment 
deposited during the drainage of Lake Agassiz [Andrews 
et al., 1995, 1999; Kenvin, 1996; Barber et al., 1999]. 
Because Lake Agassiz resided over Paleozoic carbonates, 
its waters were depleted in 14 C, with a lower 14 C/ 12 C ratio 
equivalent to an age of >310 years [Barber et al., 1999]. 
Indeed, if we apply an additional reservoir age (AR) of 
310 years to the two radiocarbon dates within the carbonate 
and sand layers to reflect the influence of dissolved Paleozoic 
carbonate on Lake Agassiz runoff 14 C/ 12 C and potentially 
enhanced sea-ice cover in the early Holocene [Jennings et al. , 
1998; Andrews et al., 1999; Barber et al., 1999; Lewis et al., 
2012], the age of this layer decreases to ^8. 3-8.0 ka, coin- 
cident with the 8.2 ka event (Figure 2a and Table SI). We use 
this revised chronology for discussion of the CT and b 1 x O sw 
records, but note that regardless of the precise timing, the 
sedimentology and color of the sediment layers tie them 
to the drainage of Lake Agassiz [Andrews et al., 1995, 1999; 
Kerwin, 1996; Barber et al., 1999; Hillaire-Marcel et al., 

2007] , 


'Auxiliary materials are available in the HTML, doi: 10.1 029/ 
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Figure 2. HU87033-017 records, (a) Age-depth relation- 
ship with blue error bars indicating 14 C dates [Andrews 
et al. , 1999]. Black is the age model discussed in the text 
with an additional 310 yr AR [Barber et al., 1999; 
Andrews et al., 1999], gray shows the age model without 
the additional AR. (b) Weight percent sand >63 pm (pur- 
ple), calcite (light green) and dolomite (light blue) 
[Andrews et al., 1999], (c) Mg/Ca calcification temperature 
(CT) with ±1.3°C uncertainty indicated, (d) N. pachyderma 
(s) /) 18 O c [Andrews et al, 1999]. (e) /) 18 O sw with propagated 
CT and d r8 O c uncertainty (±0.3%o). Horizontal dashed lines 
denote CT, <5 ls O c , and <5 18 O sw data from 49 cm core depth. 

[ 7 ] The N. pachyderma (s) CT record shows ^2°C of 
warming ~1 1.5-1 1.2 ka, with a more gradual trend to a CT 
maximum of ^8°C at ^9.7 ka (Figure 2c). CT subsequently 
cooled to 5-6°C between ~9.5 and 8.3 ka, reaching a mini- 
mum of ~3°C at ~8.3 ka. At ~8.0 ka, CT warmed to ~4°C 
and reached ~5°C by ~7.5 ka, which is equivalent to late 
Holocene CT (Figure SI and Text SI, section S7). Our 8 1 8 O sw 
record shows several deviations from the N. pachyderma (s) 
6 18 O c record of Andrews et al. [1999] (Figures 2d and 2e). 
Between ~11.2 and 9.7 ka, a 8 O c increased by ~0.4%o, 
whereas 6 18 O sw increased by ~ 1 ,2%o; thereafter, both records 
decreased by ~1.0%o. At the same core depth as the ~3°C of 
CT cooling but just preceding the increase in weight percent 
sand, <5 18 O sw further decreased by ~1.0%o. 6 r8 O sw then 
increased by ~0.9%o after ~8.0 ka, and by another ^0.4%o 
by ~7.0 ka. 

4. Discussion and Conclusions 

[8] The gradual increase in 6 18 O sw ~ 1 1 .5—9.7 ka may 
document the stabilization of the Laurentide Ice Sheet margin 
as it retreated onto the Labrador coast [Andrews et al., 1995, 
1999]. The subsequent decrease after ~9.7 ka (Figure 2b) 
likely reflects increased iceberg calving and meltwater plume 
deposition of sediment from renewed Laurentide retreat and 
the beginning of ice break-up in Hudson Strait after the 
Noble Inlet readvance, which is supported by the increase in 
weight percent sand (Figure 2b) from sediment in meltwater 


plumes and icebergs [Andrews et al, 1995, 1999; Jennings 
et al., 1998]. 

[ 9 ] The ~1.0%o further decrease in 6 18 O sw ~8. 3-8.0 ka 
identifies increased freshwater discharge into the Labrador 
Sea coincident with the drainage of Lake Agassiz. This 
depletion event was obscured in the 6 1 8 O c record by the 
contemporaneous CT cooling of ~3°C in the Labrador Sea 
that reflects the impact of cold Lake Agassiz runoff and 
regional cooling during the 8.2 ka event. The masking of this 

8 O c depletion by CT cooling is supported by climate model 
simulations of the break down in the relationship between 
temperature and 8 ls O of foraminiferal calcite during the 8.2 ka 
event [LeGrande et al., 2006], which is also observed at 
different subsurface depths in the North Atlantic with the 
application of Mg/Ca paleothermometry [Came et al., 2007; 
Thornalley et al., 2009]. Earlier studies of the Labrador Sea 
did not account for the effect of ocean cooling on 6 18 O c 
during the 8.2 ka event, explaining why a <5 18 O c decrease 
from the drainage of Lake Agassiz was previously suggested 
to be lacking [Keigwin et al., 2005; Hillaire-Marcel et al., 
2007, 2008], 

[ 10 ] The initiation of the ^8.3 ka 6 18 O sw decrease coinci- 
dent with increased weight percent calcite and dolomite 
[Andrews et al., 1999] detects the arrival of Lake Agassiz 
freshwater [Andrews et al, 1995, 1999; Kerwin, 1996; 
Barber et al., 1999; Hillaire-Marcel et al., 2007]. The 
ensuing increase in weight percent sand probably records 
iceberg rafting and meltwater plume deposition of sediment 
during the breakup of the Laurentide Ice Sheet over Hudson 
Bay (Figure 2) [Andrews et al., 1995, 1999]. Although mul- 
tiple Lake Agassiz drainage events may have occurred lead- 
ing to the final opening of Hudson Bay, the temporal spacing 
of these has yet to be quantified and was likely in a short 
interval of time [Clarke et al., 2004; Ellison et al., 2006; 
Hillaire-Marcel et al., 2007; Lajeunesse and St-Onge, 2008; 
Roy et al., 2011]. The <5 18 O sw anomaly does not show mul- 
tiple peaks that would be expected from two or more Lake 
Agassiz drainage events (Figure 3e), and thus may record the 
final drainage event or integrate several drainage events due 
to our sampling resolution relative to the hypothesized 
months-long duration of the drainage events. The 6 18 O sw 
anomaly also persisted longer than the carbonate anomaly, 
likely reflecting both the final break-up of the Laurentide Ice 
Sheet, suggested by the sand layer [Andrews et al., 1995, 
1999] (Figure 2), and sustained Lake Agassiz runoff routing 
to the Labrador Sea for several hundred years after lake 
drainage [Clark et al., 2001; Hillaire-Marcel et al., 2007; 
Carlson et al., 2009]. 

[ 11 ] We compare our <5 18 O sw decrease with other records of 
freshwater discharge during the 8.2 ka event to trace the path 
of freshwater from Hudson Bay. After flowing through 
Hudson Strait, we document the drainage of Lake Agassiz in 
the western Labrador Sea with the ~1.0%o <5 18 O sw decrease. 
In the northwest Atlantic at 43-37°N, 6 18 O c records show 
depletions of 0.4-0.6%o, reflecting wanning and/or increased 
freshwater discharge [Keigwin et al, 2005]. However, 
(> 18 O sw increases by ~0.2%o at 35°N in the northwest 
Atlantic (Figure 1) [Cleroux et al., 2012]. In contrast, (> 18 O sw 
decreases in the northeast Atlantic of ~0.8%o at ^27°W to 
~0.4%o at ~ 1 8'" W from foraminifera living at different sub- 
surface water depths (Figures 1 and 3) [Ellison et al., 2006; 
Came et al., 2007; Thornalley et al., 2009]. These records 
thus document the dispersal of Lake Agassiz freshwater 
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Figure 3. 8.2 ka event records, (a) GISP2 b 18 0 [Alley etal., 
1997], and North Atlantic <5 18 O sw records for (b) RAPiD-12- 
lKon Globorotalia inflata [Thornalley et ai, 2009], (c) ODP 
Site 984 on Neogloboquadrina pachvderma (dextral) [Came 
etal, 2007], (d) MD99-2251 on Globigerina bulloides using 
the preferred age model of the original authors [Ellison etal, 
2006], and (e) HU87033-017 on N. pachvderma (s). 


across the North Atlantic, with freshwater transported 
southwards as far as ■~37°N in the northwest Atlantic, and 
eastward into the northeast Atlantic. The decrease in b 18 O sw 
anomalies to the south and east of the Labrador Sea likely 
shows the dilution of the freshwater signal along these 
transport paths (Figure 1). 

[ 12 ] A 1 0-year simulation with the MITgcm high-resolution 
(0.167° x 0.167°) ocean-ice model suggested that a 1-year 
freshwater discharge through Hudson Strait would not reach 
deep-water formation sites and would rather be entrained into 
the subtropical gyre as far south as 25°N near the Florida 
Strait [Condron and Winsor, 2011], conflicting with avail- 
able b I8 O sw records (Figures 1 and 3) [Ellison et al., 2006; 
Came et al., 2007; Thornalley et al., 2009; Cleroux et al., 
2012]. The 10-year duration of the MITgcm simulation may 
be of insufficient length to model the transport of Lake 
Agassiz discharge from Hudson Strait to the northeast 
Atlantic. The simulations also did not include the longer- 
duration freshwater discharge from Laurentide Ice Sheet 
retreat and continental rerouting [Andrews etal., 1995, 1999; 
Clark et al., 2001; Hillaire-Marcel et al., 2007; Carlson 
et al., 2009] that could have larger effects on AMOC 
strength [Meissner and Clark, 2006; Clarke et al., 2009]. 

[ 13 ] Conversely, the UVic high-resolution (0.2° x 0.4°) 
ocean model that includes a simplified atmosphere found that 
freshwater on the Labrador shelf could affect the AMOC 
within years of its discharge [Spence et al., 2008]. In addi- 
tion, the UVic model simulated that the duration and maxi- 
mum amplitude of the freshwater forcing and AMOC 
response are relatively insensitive to increasing resolution 
from low-resolution general circulation models (GCM) to 
high, eddy-resolving resolution [Spence et al., 2008], 
We also find agreement between changes in b 18 O sw by the 
NASA Goddard Institute for Space Studies fully-coupled 


GCM ModelE-R that includes water isotopes throughout the 
hydrologic cycle [LeGrande and Schmidt, 2008] and the 
observed decrease in b 18 O sw records during the 8.2 ka event 
(Figure 1), similar to other coupled climate model studies 
[Meissner and Clark, 2006; Spence et al. , 2008; Clarke et al. , 
2009]. These model-b 18 O sw comparisons suggest that 
although Lake Agassiz runoff may have been initially 
trapped in boundary currents along eastern North America 
[Keigwin et al., 2005; Wunsch, 2010; Condron and Winsor, 
2011], the freshwater eventually escaped from the continen- 
tal shelf, was entrained in the Gulf Stream and North Atlantic 
Current, and reached regions of deep-water formation in the 
northeast Atlantic where it could affect the AMOC [Ellison 
et al., 2006; Kleiven et al., 2008; Hoogakker et ai, 201 1], 

[ 14 ] In conclusion, we document substantial near-surface 
cooling of the Labrador Sea during the 8.2 ka event. After 
accounting for this cooling, our <5 18 O sw record shows a sig- 
nificant decrease due to the drainage of Lake Agassiz; 
a signal that was previously masked in b 18 O c records by the 
near-surface cooling. When combined with b 18 O sw records 
from elsewhere in the North Atlantic, we can trace the path- 
way of Lake Agassiz discharge from Hudson Strait to regions 
of deep-water formation in the northeast Atlantic, supporting 
the hypothesis that the drainage of Lake Agassiz and atten- 
dant runoff routing into the Labrador Sea forced the 8.2 ka 
event [Klitgaard-Kristensen et al., 1998; Barber et al., 1999; 
Clark et al. , 200 1 ; Teller et al. , 2002; LeGrande et al. , 2006; 
Hillaire-Marcel et al., 2007; Clarke et al., 2009; Carlson 
et al., 2009]. 
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